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THE LEADING- EDGE EFFECTS ON THE 
LAMINAR FLAT-PLATE BOUNDARY LAYER AND THE 
AERODYNAMIC HEATING AT MACH 10.4* 

By Howard W. Stone 
Langley Research Center 

SUMMARY 

Wall-pressure and heat-transfer distributions and boundary -layer impact-pressure 
profiles were obtained on a flat plate with both sharp and cylindrically blunted leading 
edges. The range of free- stream Reynolds number based on leading- edge thickness was 
from 73 to 70 000. Although the boundary-layer displacement effects behind a sharp 
leading edge are known to increase the wall pressure and heat transfer in the weak- 
interaction regime, these tests show that boundary- layer displacement effects result in 
only a slight change in the impact-pressure and velocity profiles. Slight blunting of the 
leading edge had little effect on the pressure distribution but, with the boundary layer 
growing in an entropy gradient, the heat transfer was reduced. Large leading- edge blunt- 
ness increased the heat transfer relative to the sharp leading edge. 

INTRODUCTION 

The effect of leading- edge thickness on the viscous flow field adjacent to a body in 
a hypersonic stream has been of interest to many researchers who have, in general, been 
concerned with the extremes of the problem; that is, the sharp leading edge or the very 
blunt leading edge. Experimental and theoretical studies of the sharp- leading- edge flat 
plate in a hypersonic flow with a laminar boundary layer are presented in references 1 
to 4. The viscous-induced weak- interaction regime is discussed in references 1 to 3, 
whereas reference 4 reviews much of the work in the very near leading-edge regime. 
Experimental and theoretical studies for both the sharp-leading-edge and blunt-leading- 
edge flat plates with a laminar boundary layer are found in references 5 to 13. Generally, 
the wall-pressure distributions are reported to be in fair agreement with the various 

*The material presented herein is based in part upon a thesis entitled "The 
Leading-Edge Effects on a Laminar Flat- Plate Boundary Layer at Mach 10.4" submitted 
in partial fulfillment of the requirements for the degree Master of Aerospace Engi- 
neering, University of Virginia, Charlottesville, Virginia, May 1968. 



weak- interaction theories when the viscous-induced effects are predominant (sharp 
leading edge) and in fair agreement with the modified blast-wave theory when the inviscid 
effects are predominant (very blunt leading edge). 

Measured laminar heat- transfer distributions for sharp- and blunt-leading- edge 
flat plates are presented in references 8 to 11. Bluntness generally tends to increase 
the heat transfer, but in reference 10 a slight decrease was shown prior to the increasing 
trend as the leading edge was progressively blunted. 

Creager (ref. 14) conducted impact-pressure surveys of the region between the 
wall and the shock wave on flat plates with various leading-edge thicknesses at Mach 5.7. 
No attempt, however, was made to compare the surveys with theory. 

The purpose of this investigation was to examine the leading-edge effect on the 
boundary -layer profiles and the boundary- layer thickness in addition to the wall pres- 
sures and heat transfer. These measurements were obtained on a flat plate behind a 
sharp leading edge and two blunt leading edges at a free- stream Mach number of 10.4. 

The impact-pressure profiles have been converted to velocity profiles and are compared 
with a zero-pressure-gradient boundary-layer solution and a boundary-layer theory using 
the local similarity concept to account for the effect of the pressure gradient. The tests 
were conducted in air with a range of free-stream Reynolds number based on leading- 
edge thickness from 73 to 70 000 for leading-edge thicknesses of 0.005, 0.16, and 
1.27 cm. 


SYMBOLS 


a speed of sound 

Cm heat capacity of the metal 

d leading-edge diameter 

h enthalpy 

K 3 pressure gradient correction factor (see ref. 12 ) 

M Mach number 

Nst^oo Stanton number based on free-stream conditions 
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p 


pressure 


2 impact pressure 

unit Reynolds number based on free- stream conditions 

Rqo d Reynolds number based on free-stream conditions and leadii^-edge diameter 

^oo Y Reynolds number based on free-stream conditions and distance from leading 

edge 

T temperature 

u velocity 

u' effective velocity 

X distance from leading edge of sharp model 

y distance perpendicular to model surface 

z lateral distance on model surface measured from model center line 


y ratio of specific heats 

6 boundary- layer thickness 

6' boundary- layer characteristic thickness 

Voo Blasius similarity variable, x 

fi viscosity 



Subscripts: 

e boundary -layer -edge condition 
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t total (stagnation) condition 

w wall condition 

oo free-stream condition 

APPARATUS AND TESTS 
Facility and Test Conditions 

The present tests were conducted in the Mach 10 test section of the Langley 
continuous-flow hypersonic tunnel. In this facility the prescribed test conditions are 
continuously maintained in the 79-cm-square test section by recirculating the test air 
through a series of compressors and electrical resistance tube heaters which heat the 
test air to avoid liquefaction. A more complete description of the facility is available in 
reference 15. 

The tunnel Mach number distributions obtained from impact -pressure surveys are 
presented in figure 1 for the nominal test stagnation chamber pressures of this study. 
The nominal test conditions obtained by assuming an isentropic expansion (with correc- 
tions for intermolecular force effects, ref. 16) are summarized in the following table: 


kN/m2 


2410 

5170 

8270 


The maximum deviation from the nominal test stagnation-pressure values for dif- 
ferent runs was 0.8 percent whereas the test stagnation temperatures varied from 966° K 
to 1040° K with no more than an 8° variation during a given run. These deviations in the 
stagnation conditions along with slight tvinnel Mach number gradients result in a Mach 
number deviation of ±0.10. 



Free-stream unit 

M 

Reynolds number 

■‘’■‘oo 

per cm, 


RcoA 

10.3 

1.5 X 104 

10.4 

3.1 

10.4 

5.1 


Models and Instrumentation 

The flat-plate models consisted of interchangeable instrumented plates supported 
by a frame with a 20° bevel extending across the leading edge and along both sides of the 
plate. (See fig. 2.) The pressure plate was 0.476-cm-thick stainless steel with 0.178-cm 
pressure orifices located as shown in figure 3. The wall-pressure measurements at 
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36 orifices were obtained by three ionization gages connected through three 12-port 
scanning valves. All the data were recorded automatically on magnetic tape by an 
analog-to-digital converter. Random instrument errors plus recording errors result in 
an accuracy of about ±6 percent of reading. 

The measured mean thickness of the inconel heat-transfer plate was 0.081 cm 
±5 percent. The material density was 8.41 x 10^ kg/m3 and from a curve fit of values 
published in the literature, the heat capacity in J/g-°K varied as 


C 


m 


0.0000611 




+ 0.0980 


4.187 


( 1 ) 


The temperature on the back side of the thin inconel plate was measured by 30-gage 
chromel-alumel thermocouples spotwelded to the plate at the locations shown in figure 4. 

The leading edges of the sharp plates were 0.005 ± 0.002 cm thick and were modi- 
fied by adding cylindrical leading edges 0,16 cm and 1,27 cm in diameter for the blunt- 
model tests. Both the heat-transfer plate and the pressure plates had mirror finishes 
with a maximum waviness of 0.33° on the center line. Detachable end plates with sharp 
leading edges were constructed from 0.32-cm-thick stainless steel. The shape of the 
end plate is shown in figure 4, 


Boundary- Layer Probes 

Two boundary-layer probes were constructed for these tests. Each probe had 
three tubes formed from stainless-steel tubing having 0,23-cm i.d. and 0.32-cm o.d. 
which was flattened and filed at the tip to form almost rectangular orifices. The three 
tubes were mounted on a 1.27-cm o.d. stainless tube which was air-cooled, A sketch of 
the probe heads and the measurements of each orifice are presented in figure 5. A 
photograph of one of the probes is shown as figure 6, 

The tube heights were less than 15 percent of the expected minimum boundary-layer 
thickness which is consistent with criteria obtained from previous probe work (e.g., 
refs. 2 and 17). Also, the tubes were flattened to a width-height ratio of approximately 7 
after the "viscous effect" study of MacMillan (ref, 18), 

Both probes were electrically driven from outside the tunnel. During bench tests 
of the traversing mechanisms, the probe location always repeated to within 0,008 cm when 
within 2.5 cm of the model surface. The surveys were conducted at the locations shown 
in figure 3. Each tube of the boundary- layer probes was connected to multiple trans- 
ducers to provide an accuracy of ±3 percent over the range of impact pressures incurred. 
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Tests 


In all tests, the model was injected into the tunnel in less than 1 second after the 
proper test conditions had been established, and the desired data were measured. Upon 
completion of these measurements, the model was withdrawn from the stream and cooled 
by air jets to nearly room temperature. Oil-flow tests were conducted on the sharp- 
leading-edge model initially to establish which end-plate configuration would provide 
parallel flow on the model center line. In the wall-pressure tests and the boundary- 
layer surveys, the pressures were monitored continuously and were recorded on mag- 
netic tape after they had reached steady state. During the surveys, when the probe 
touched the model surface the fouling circuit signaled and a data point was recorded to 
establish the surface location on the recorded data. The probe was moved only in the 
direction toward the model surface throughout all surveys. In the heat-transfer tests, 
the temperature at each thermocouple was recorded 20 times per second for about 
5 seconds. The Stanton number was obtained by applying the transient- calorimeter tech- 
nique which involves calculating the rate of heat storage in the thin metal skin from the 
slope of the measured temperature-time data as described in reference 15. In the pres- 
ent tests, free- stream flow properties were used throughout the calculation. 

DISCUSSION OF RESULTS 

In this investigation, an end plate was installed on the strut side of the model and 
extended upward at an angle of 30^ to the plane of the model surface. In this configura- 
tion, the model surface was shielded from the interference field of the strut. The oil- 
flow pattern (see fig. 7) shows that except for a slight inward flow near the leading edge 
on the far side and some outward flow along the sides, the flow over the plate was parallel 
to the center line. The spanwise pressures for this configuration at the 90-percent- chord 
location were constant within the accuracy of the data. 

Pressure Distributions 

The bluntness and viscous-interaction effects on the wall-pressure distributions 
are shown in figure 8 where along the model center line is plotted as a function 

of X for three unit Reynolds numbers. The modified blast-wave theory of Bertram 
and Blackstock (ref. 13) is generally in good agreement with the pressures behind the 
1.27-cm-diameter leading edge where the bluntness- induced effects are dominant. There 
is, however, some difference between data and theory on the rear of the model due prob- 
ably to transducer inaccuracy. The data for the models with the 0.16-cm-diameter and 
the sharp leading edges agree almost within the accuracy of the data. 
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The wall pressures for the models with the sharp- and 0.16-cm-diameter leading 
edges tend to increase slightly on the rear third of the model but the oil-flow patterns 
and the lateral pressure distributions did not indicate any model edge effects in this 
region. There does appear to be a model trailing- edge effect over the last 5 cm of the 
model. 

The "complete" theory of Bertram and Blackstock (iteration of eqs. (7) and (15) of 
ref. 13) for predicting viscous -induced pressures, shown in figure 8, does not predict the 
sharp- leading- edge data. The theory assumes an isothermal wall ^T-^/Tj; « 0.34 in the 
present tests^ but the temperature of the leading edge itself will approach the stream 
total temperature which may have an effect on boundary- layer growth and thus on the 
induced pressure rise. Also a slight misalinement of the model with respect to the tun- 
nel flow may account for most of this difference at the highest unit Reynolds number; 
however, the greater part of the difference at the lower unit Reynolds number is 
unexplained. 


Impact- Pressure Profiles 

Changes in the boundary layer and the adjacent flow field due to changing the 
leading- edge thickness are shown in the measured impact-pressure profiles presented in 
figure 9. Comparison of the impact-pressure profiles for the models with the sharp- and 
the 0.16-cm-diameter leading edges reveal a marked leading-edge effect whereas the 
wall-pressure distributions were very similar. The profiles at x = 71.6 cm for the 
1.27-cm-diameter leading-edge model (fig. 9(c)) were reduced by using the blast-wave 
values of p^ instead of the measured values. The profiles at both probe locations for 
the 1.27-cm-diameter leading-edge model agree outside of the boundary layer where the 
impact pressures were equal to the values calculated by assuming an isentropic expan- 
sion from the model stagnation line. Thus, the boundary layer behind the 1.27-cm- 
diameter leading edge grew in a region of constant entropy associated with the normal 
portion of the bow shock wave. 


Velocity Profiles 

The velocity profiles for the sharp-leading- edge model at each unit Reynolds num- 
ber are presented in figure 10. These data were obtained from the measured impact- 
pressure profiles by using a theoretical total- enthalpy distribution due to the theory of 
Klunker and McLean from reference 19 (discussed in the appendix). As the figure shows, 
the near-wall data points are affected by the wall-probe interaction and possibly by the 
low Reynolds number effect. These data points are neglected, therefore, throughout the 
rest of the discussion. The theory of Klunker and McLean and the data are in good agree- 
ment at all Reynolds numbers. Thus, a zero-pressure-gradient boundary-layer theory 
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based on free- stream conditions will provide, for these test conditions, a satisfactory 
estimate of the boundary-layer characteristics for 0.57 = Xoo = which is assumed 

in the weak-interaction theory. 

The wall pressures on the sharp-leading- edge model are in a range where the mea- 
surements presented in reference 20 show that there is a small effect of low pressure on 
the viscosity. Therefore, the Klunker and McLean theory was also calculated using the 
Kurdt and Warberg equation (ref. 20) for the variation of viscosity with pressure. The 
resulting theoretical velocity profiles were only slightly different from the theoretical 
velocity profiles presented in figure 10 . 

The velocity profiles for the 0.16-cm-diameter leading-edge model, taken at 94 and 
176 diameters from the leading edge, are presented in figure 11 along with the profiles of 
the Klunker and McLean theory for a sharp leading edge for comparison. The theory 
also provided the total- enthalpy distribution that was used to obtain the velocity from the 
measured data. The data show that there is a small velocity gradient outside of the 
boundary layer. In figure 11(a), for x = 38.1 cm, the "knee" of the profile normally asso- 
ciated with the boundary-layer edge appears to be at 77 ^ = 32. This location corresponds 
to y = 1.53 cm on the impact-pressure profiles presented in figure 9(b) or ^p^ 2 /^w) 

is approximately 27 at the boundary- layer edge. The impact-pressure gradient outside 
of the boundary layer indicates that the boundary layer is growing in an entropy gradient. 

Behind the 1.27 -cm- diameter leading edge, the boundary layer has been shown to 
be growing in the region of constant entropy adjacent to the model surface. The local 
similarity theory discussed in the appendix was used to calculate velocity profiles for 
this configuration by assuming that the Mach 10.4 normal-shock total pressure was the 
boundary-layer-edge total pressure. The velocity profiles were computed from the mea- 
sured impact-pressure data by using the total -enthalpy distribution calculated by the local 
similarity theory. The theory and the velocity-profile data shown in figure 12 are in good 
agreement at x = 38.1 cm but do not agree as well at x = 71.6 cm. 

To obtain an indication of the bluntness effect on the velocity profiles, a boundary- 
layer characteristic thickness 6 * was defined as the point where a straight-line fairing 
of the linear portion of the boundary-layer velocity profile crosses a straight-line fairing 
of the velocity data points outside of the boundary layer. This point also defines an effec- 
tive boundary-layer edge velocity Ug. A plot of u/ug as a function of y/ 6 ' for the 
three leading-edge configurations is presented in figure 13. The bluntness effect is con- 
fined to the outer portion of the boundary layer where the boundary layer is thickened and 
the velocity gradient decreased as the leading-edge thickness is increased. 
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Heat Transfer 


The heat-transfer data, presented in figure 14, show that progressively blunting 
the leading edge first decreases and then increases the wall heat transfer relative to that 
for the sharp-leading- edge model. Neal (ref. 10) found a similar bluntness effect at 
M^ = 6.8. 

The results of the Bertram and Feller (ref. 12) method to "correct" the heat- 
transfer data for the effect of the induced pressures are shown in figure 15. Both the 
"corrected" and "uncorrected" heat-transfer data for the sharp-leading-edge model are 
compared with the zero-pressure-gradient theory of Klunker and McLean. The fairing 
through the uncorrected data is parallel to and 19 percent above the theory. There is 
good agreement between the corrected data and the theory. 


Boundary- Layer Thickness 

Creager (ref. 21) found fair agreement between his Moo = 4.0 boundary-layer- 
thickness data and the equation 



1.73 /Tw\ 
Moo^VtJ 


+ 0.332(y - 1) + 



( 2 ) 


Comparison of the present test results with this equation required locating the boundary- 
layer edge, usually defined as the location of u/ue = 0.999. However, this definition is 
difficult to apply to the intermediate bluntness profiles with the velocity gradient outside 
the boundary layer. In an attempt to select a comparable boundary- layer edge value, a 
straight-line fairing of the velocity data outside the boundary layer was assumed to be 
analogous to the u/ug = 1.0 line. The tangent point of the data with this fairing was 
selected as the boundary-layer edge. These data are compared with a plot of equa- 
tion (2) in figure 16. The Monaghan T' prediction (ref. 22) based on the free-stream con- 
ditions is also shown in the figure. Both methods provide a satisfactory estimate of the 
boundary-layer thickness for the sharp- and 1.27-cm-diameter leading-edge data and a 
portion of the 0.16-cm-diameter leading-edge data. Creager's data (ref. 21) showed that 

the boundary -layer thickness is proportional to the square root of x for 0.6 = ^ S 15.5 

d 

at Mach 4. The present data show that 5 <x \Jx for 30 i ^ 14 000 at Mach 10.4. 

d 


CONCLUSIONS 


The investigation of the leading- edge effects on the laminar boundary layer of a 
flat plate with a sharp and two cylindrically blunted leadir^ edges in a Mach 10.4 air- 
stream support the following conclusions: 
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1. The boundary-layer solution presented in NACA TN 2916, based on the assump- 
tion of a constant wall pressure, was in good agreement with the velocity-profile data 
for the sharp flat plate with a weak- interaction wall -pressure distribution. The heat- 
transfer data agreed with the theory after a modification for the wall pressures by the 
method presented in NASA MEMO 5- 24- 59 L. 

2. The profile data for the 0.16-cm-diameter leading-edge model, taken at 94 and 
176 diameters from the leading edge, indicated that a large impact-pressure gradient 
and a small velocity gradient existed outside the boundary; that is, the boundary layer 
was growing in an entropy gradient. This bluntness thickened the boundary layer and 
slightly reduced the heat-transfer coefficients relative to the sharp- leading- edge data. 

3. The velocity-profile data for the 1.27-cm-diameter leading-edge model were in 
fair agreement with a local similarity boundary-layer theory. The large bluntness 
increased the heat-transfer coefficients relative to the data for the sharp- leading- edge 
model. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 14, 1969, 

129-01-07-07-23. 
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APPENDIX 


ANALYTICAL METHODS 
Laminar Boundary -Layer Theories 

The theory of Klunker and McLean (ref. 19) is an easily programed iterative solu- 
tion for the two-dimensional compressible laminar boundary -layer equations in the 
absence of a pressure gradient and for an isothermal wall. The solution includes the 
temperature variation of all fluid thermal properties, and the results are representative 
of the constant-pressure solutions for a flat plate. The Blasius similarity parameter is 
used to reduce the boundary -layer equations for steady flow along an isothermal wall to 
total differential equations. After the equations are nondimensionalized and simplified 
by integrating factors, a rapidly converging, successive approximation, integration 
method is used to obtain a solution. In the present application the fluid thermal prop- 
erties are obtained from the NBS tables (ref. 23). The work of Grieser and Goldthwaite 
(ref. 20) is used to evaluate the viscosity. The equation for viscosity in N-sec/m2 

t3/2 ft 

M = 1.419 ^ 10-6 (Al) 

T + 99.22 

represents the best fit of their data over the temperature range from 80° to 294° R. In 
the present investigation the free-stream conditions were used throughout the solution. 

The local similarity concept can be used to account for the effect of the pressure 
gradient on the boimdary layer. A programed numerical solution for the locally similar 
equations of reference 24 was used to obtain theoretical estimates in this investigation. 
The solution used was for an ideal gas with a Sutherland law form of viscosity variation. 
This solution requires prior knowledge of the boimdary -layer edge flow properties and, 
therefore, can be used satisfactorily only when the static- and total -pres sure distributions 
are known over the entire length of the model. The wall was assumed to be isothermal. 

Total -Enthalpy Profiles 

To obtain velocity profiles from the measured impact -pressure profiles, it is nec- 
essary to know the corresponding total -temperature or total -enthalpy profile. Unfortu- 
nately, there was no reliable total -temperature probe available that could be used in 
these tests at the time the tests were performed and, therefore, an appropriate theoret- 
ical profile was used. Any error due to this assumption is reduced in the velocity pro- 
files, however, since the velocity is proportional to the square root of the difference 
between the total and static enthalpy and the values are fixed where the difference 
approaches zero. 
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APPENDIX 


The calculation procedure involved applying the energy equation 

ht = h + (A2) 

to each data point in the boundary -layer profile. At each point M was known from the 
measured profile and an appropriate value of h^ was known from the assumed total - 
enthalpy distribution. Equation (A2) can be solved either by assuming an ideal gas and 
using the appropriate definitions of h and a or by an interation procedure using tables 
of values of h and a as a function of temperature such as the NBS tables (ref. 23). 
Both procedures were used and yield almost the same results for the conditions of the 
present tests. 
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Model surface location 


Figure L- Tunnel Mach number calibration. 
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Figure 2.- Model positioned in test section. L-65-5823 
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Figure 3.- Pressure-orifice and boundary- layer probe survey locations. Ail dimensions are in centimeters. 
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Figure 6.- Photographs of boundary-layer probe head. 
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Figure 7.- Oil-flow patterns. 
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Figure 10.- Velocity profiles for sharp- leading-edge model. 
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Figure 16.- Boundary- layer thickness behind different leading edges. 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 


0 2 U 0 J 1. 3 7 
Air; FjiiCr- 

Kl'<VLAr.:r; r 


5 1 30 S 

capons ladopai o 
R iORCc TASh» u 


Rnoo 0090 3 
a / AFWL/ 

A'tXloO 87 11? 


A T 1 


e . 


0 0 


A c r i ■ 


r T I 


L 1' 


POSTMASTER; 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


”Tbe aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof/^ 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 

technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published i'n a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washingion, D.C. 20546 



